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Abstract 

The canonical transformation diagonalizing the one-particle tight binding 
Hamiltonian for an alternating chain with two non-equivalent sites per unit 
cell has been used to introduce the Hubbard-like interactions (on-site, inter- 
site, bond-site and intra-bond) in the corresponding two-band model. The 
considerations have been restricted to the particular case, suggested by the 
undistorted 3/4-filled CuOj, chain, where the alternating structure comes only 
from non-equal on-site atomic energies and the gap between the two bands is 
sufficiently large. The renormalization group method has been then applied 
to the upper band of this model and the corresponding phase diagram has 
been analyzed in terms of some renormalized (density dependent) Hubbard- 
type couplings for arbitrary filling of the upper band. The particular case of 
CuO% chain is also discussed and the corresponding phase diagram has been 

drawn in terms of the original (and realistic) coupling constants. 
PACS numbers: 71.10 Fd, 71.10 Hf, 71.10 Pm, 74.20 Mn 
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I. INTRODUCTION 



A large amount of work has been devoted to the study of alternating chains not only 
due to their interesting structure, but also in connection with a better understanding of the 
high-T c superconductivity mechanism The same phenomenon renewed the interest in 
the Hubbard model and, quite recently, in its one-dimensional extensions |1. However, 
the particularities induced by the Hubbard-like couplings in alternating chains are less well 
known: the charge gap that opens in a one-dimensional dimerized Hubbard model has 
been estimated in various limits using analytical results and exact diagonalizations of small 
clusters |§; a model where the dimerization is induced via alternating on-site repulsions has 
been studied within a boson representation theory ||] and also by renormalization group 
(RG) technique f5|. The aim of this paper is to investigate, using the RG method, the effect 
of the Hubbard-type interactions on the ground-state properties of a chain with alternating 
on-site atomic energies. 

The paper is structured as follows. The one-particle Hamiltonian in the tight binding 
approximation corresponding to an alternating chain with two nonequivalent sites per unit 
cell can be diagonalized by a canonical transformation; one gets a two band model. The 
Hubbard-type interactions, i.e. interactions which in the site representation couple only 
electrons belonging to the nearest neighbor sites, give rise to both intra- and inter-band 
couplings; however, if the gap between the two bands is sufficiently large and the system 
is more than half-filled, as for the CuO^ chain occurring in high-T c superconductors, the 
last ones can be neglected in describing the low energy physics. The obtained expressions 
of the Hubbard-type interactions (upper band), in the particular case of alternating on-site 
energies and equal hopping amplitudes, close Sec. II. The standard RG analysis (second 
order) ||, briefly reviewed in Sec III, is done in terms of the g-constants describing the 
elementary processes of forward, backward and umklapp scatterings: their expressions are 
obtained by evaluating the Hubbard-type interactions (upper band) at the Fermi points. 
Using the scaling to the exact soluble models Tomonaga-Luttinger (TL) |7|] and Luther- 
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Emery (LE) ||, we can predict the low energy physics of our system; the ground-state 
phase diagrams in terms of the model parameters and at arbitrary band filling are presented 
in Sec. IV. The relevance of our analysis to the case of the undistorted 3/4-filled CuO^ chain 
and the corresponding phase diagram in terms of the original coupling constants (neglecting 
the intra-bond coupling which usually is very small) are discussed in Sec. 5. The last section 
summarizes the main results of this work. 



II. THE TWO-BAND MODEL 

Let us consider the alternating chain from Fig. 1, with two non-equivalent sites per 
unit cell (N cells, periodic boundary conditions). Within the tight binding approximation, 
the second quantized form of the one-particle Hamiltonian in the site representation of the 
atomic orbitals reads 101 
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-E[(ta} + *a} +1 )& i + J7.c.] (1) 



where j = 1,N is the cell index, and the a(6)-operator corresponds to the annihilation of 
one electron on a site A(B); we ignore for the moment the spin variable a =|, \ of the 
electron (it can also be considered as included in the cell index). The constants e and t 
(t) are positive and denote respectively the site energies and the parameter describing the 
hopping between A and B sites belonging to the same cell (nearest neighbor cells); all of 
them can be expressed in terms of the atomic quantities ||. From the site representation 
we can pass to the momentum representation by a usual Fourier transform; Eq. (1) becomes 
then 



^o = eE ( a l a k ~ b lh 

k 

-Y,[{t + te- iak )a{b k + H.c] (2) 



where k G (— ir/a, ir/a] with a being the lattice constant. 



By mixing the a- and 6-operators, the Hamiltonian H,q has not a diagonal form; it can 
be diagonalized by the canonical transformation [Bf 



with 



A(k)e-^c hk + B(k)c 2 , 



b k = B(k)c lk - A{k)e i4, ^a 



2. A- 



A(k) 



V2 



e(k) 



m = 7=2 



1 + 



e(k) 



tan[0(&)] 



e(k) in Eq. (4) has the form 



t sin(afc) 
t + tcos(ak) 



e{k) = J A 2 +4tt cos 2 (ak/ 2) 



(3) 



(4) 



(5) 



where 



A = Ve 2 + (t - I) 2 
In terms of the c-operators, the one-particle Hamiltonian (2) reads as 

n = J2(-l) a e(k)c{ )a c kja , a = 1,2 

A. a 



(6) 



(7) 



and defines the kinetic part of the two-band model: the dispersion law in the upper (lower) 
band is given by plus (minus) e{k) with a gap between the two bands equal to 2A; in the 
limit A — » , one recovers the usual dispersion law —It cos(ak) for an ideal (non-alternating) 
structure. 

By analogy with the usual case, the Hubbard-type interactions between the electrons in 
an alternating chain are introduced in the site representation as follows: 



(i) on-site 



3,0 



where nf}f ] = a\ a a^{b\ a b^)- 
inter-site 



3W 



(iii) bond-site 

^3 = E Ha { X Ah„ + ^-1 J + H - C -_ 

+ E (x B a jia + X B a j+lia ) + E.c. 

(iv) exchange-hopping 

^ = -\ E h> (y&iA. 



n 



n 



3 -a 
B 



3 -a 



(v) pair-hopping 



ns = lY,HA-«( wb 3,-A<, 



(8) 



(9) 



(10) 



(11) 



(12) 



The expressions of the Hubbard-type interactions in the corresponding two-band model 
are obtained by passing in Eqs. (8)-(12) to the momentum representation and replacing 
after that the a k - and 6 fc -operators by the c a fc -operators, according to Eq. (3). Below we 
shall restrict our considerations only to the particular case of alternating on-site energies 
(e 7^ 0) and equal hopping amplitudes (t — t); in this case (ai = a 2 = a/2) 



ak 



(13) 



and there will be no distinction between intra-cell interaction constants (without overline) 
and the inter-cell ones (with overline). Any interaction in terms of the a- and 6-operators 



gives rise to both intra- and inter-band couplings. Nevertheless, if the gap between the two 
bands is large enough (comparable to the bandwidth and larger than all electron couplings) 
the considerations can be restricted to the partial filled band; assuming this fact together 
with a concentration greater than one electron per site, we shall consider further only pro- 
cesses from the upper band (the realistic case of the Cu0 3 chain will be discussed in Sec. 
V). Our (one-band) model Hamiltonian reads thus as (the index a = 2 will be omitted in 
c-operators) 

t 1 5 

^ = ^ £ (k) C k,a C k,a + T^y X! X! $k 1 +k2,k:i+k4 
k,a (5=1 ki-4,;a,a' 

x Vpih, .., h; a, <j')cl ua cl 2 a ,c k4>a ,c k3ta (14) 



where e(k) is given by Eq. (5); Vg quantities ((3 = 1,5) correspond, respectively, to the five 
types of the Hubbard interactions defined by Eqs. (8)-(12), and they have the following 
expressions: 



Vi = 



U A \{B{k l )±{A^B) 



i=l 



S a ,-a' (15) 



V 2 = 2Vcos[a(h - fc 3 )/2] 

x [Aik^B^Aik^Bih) ± (A <-> B)} (16) 

V 3 = -A{X A [cos(ak 1 /2)A(k 1 )B(k 3 ) + cos(aA;3/2) 

x B{k l )A{h)] B(k 2 )B(k A ) ±(A^B)} 5^ a , (17) 

V 4 = 2Wcos[a(k 1 - fc 4 )/2] 

x [AikjB^Bik^Aih) ± (A <-> B)} (18) 

V 5 = 2Wcos[a(k 1 + k 2 )/2] 

x [A{k 1 )A{k 2 )B{k^)B{k i ) ±(A<-> B)} 5 a ,- a , (19) 

The 5-function in Eq. (14) assures the conservation of the total momentum up to a reciprocal 
lattice vector, i.e. k\ + k 2 = k 3 + k 4 + Q with Q = or Q = ±2n/a. In Eqs. (15)-(19), the 
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upper (lower) sign corresponds to the normal (umklapp) scattering, i.e. to Q = (±27r/a); 
this fact, characteristic of an alternating structure ||, comes from the phase factor </> of the 
canonical transformation (3). Let us also note that for the lower band (a = 1), V2, V4 and 
V5 have the same expressions while V\ — > V\ [Ua Ub] and V3 — > — V3 [Xa <-> Xb\- 

III. RENORMALIZATION GROUP ANALYSIS 

The low energy physics of our model (14) can be described, within the RG method ||, 
by assuming: (i) a linear dispersion law 

e(jfe) ~ v F (k F - \k\) , v F > (20) 

together with the existence of a momentum cut-off k restricting all possible states of the 
electrons to those around the Fermi points (or equivalently, a bandwidth cut-off E = 2t>^/c ); 
(ii) all interaction processes can be classified into four different types with the coupling 
constants g i (i = 1,4) obtained by evaluating the bare potentials at the corresponding 
values + or —hp of the momenta. 

With the usual notations 0, the expressions of the (^-constants corresponding to the 
model (14) and also their form for a non- alternating chain are given in Table I, where we 
have introduced the renormalized Hubbard constants 

U = U A B\k F ) + U B A\k F ) 
V = 2VA 2 {k F )B 2 {k F ) 

(21) 

X = X A A(k F )B 3 (k F ) + X B A 3 (k F )B(k F ) 

W = 2WA 2 {k F )B 2 {k F ) 

Unlike for a non- alternating chain, they now depend on the band filling n [ak F = 7r(l — n), 
< n < 1]. As can be remarked from Table I, there is an almost complete analogy between 
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the case of an n- filled band coming from an alternating chain and a (n+l)/2-filled band of a 
non-alternating one; the differences consist in the renormalization of the Hubbard constants 
and, obviously, in the umklapp process. 

All physical results predicted by the RG method (second order) can be discussed in terms 
of four independent coupling constants: g 1± , g 3 , 



9a = 0111 - 0211 + 92± 



(22) 



and 



9 P = 0111 - 0211 - 9 2 ± 



(23) 



(the main effect of g A can be included in a renormalized Fermi velocity). The coupling 
constants g a and g 1± describe the spin sector; g p and g 3 , the charge sector. The RG equations 
(in units of ttvf) read as || 

r d 9., 



dx 



(x)=g 2 1± {x)[l + lg a {x)} 



d 9i± 
- dx 



(*) = 9*&)g 1± (x)[l + \g a {x)\ + \cftjx) 



(24) 



f ^)=9l(x)[l + lg p (x)] 



= 9 p (x)g 3 (x)[l + \g(x)\ + \gl{x) 



(25) 



k dx r, ^v-/^v-/i- ■ 4^pv 

and reflects the charge-spin separation; x = ln(E/E ) G (— oo,0], with E a smaller cut-off 
than the original one E . By solving Eqs. fl2~4| ) and (|25]), a set of equivalent problems, 
related by RG transformations, can be found; the low energy physics is essentially the same 
for all models with the g-constants along a certain g(x) solution. The corresponding flow 
diagrams indicate two distinct regimes: (i) for g a > \g 1A _\ (g„ > \g 3 \) the spin (charge) part 
of the system scales to the exact soluble TL model [|7| which is also the fixed point (this is 
the weak coupling regime); (ii) for g a < \g 1± \ (g„ < \g 3 \) the spin (charge) part of the system 



S 



scales to one of the two strong coupling fixed points, but before this it crosses the LE line 
and thus its behavior can be inferred from the exact solution || . 

Let us note that for the Hubbard-type interactions, due to the S77(2)-spin symmetry of 
the Hamiltonian (14), g a is always equal to g 1J _; consequently, we get four distinct regions 



/ : < 



9a >0 

9 P > \g 3 \ 



II : < 



9„>0 
9 P < \9 Z \ 



(26) 



9a < I 9a < 

III : { IV: 



9n < Ifl-sl 



9 P > \g 3 \ 



where, in our case, 



g a = U - 2V cos(t™) - 8X sin(7m/2) + AW 



(27) 



g p = -U - 2V[2 + cos(Trn)] + 8X sin(7rn/2) 

+4Wcos(t™) (28) 

and g 3 can be read from Table I. In terms of the model parameters, the four regions defined 
by Eq. (26) can be conveniently described in the (W, V)-plane as follows: 

W > F | W > F 

1 :{ ~ II: ~ 

V < G V > G 



(29) 



W < F \ W < F 

I I I : { IV: I 

V > G V < G 



where 



F — \—U + 2V cos(vm) + 8X sin(7rn/2)] /4 (30) 
does not depend on W and 
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f 8;r sin(7m/2) + 4Wcos(vm) - U 



2 [2 + cos(vm)] 



G 



(31) 



^Max{Y 1 ,Y 2 } 



is independent of V where 



Y x = U A B 4 - AV2X A AB 3 , Y 2 = Y 1 [A^ B] 



(32) 



A and B in Eq. (32) stand respectively for A{k F ) and B{kp) [see Eq. (4)] evaluated at 
half-filling. The scaling of our system to one of the two exactly soluble models, TL or LE, in 
each of the four regions, at any density and for both charge and spin sector is summarized 
in Table II. 



Based on the exact results obtained for the TL and LE models and using the scaling 
arguments, we can now describe the low energy physics of our system. Anytime the system 
scales to the LE model, there is a gap in the corresponding charge or spin sector; in the 
TL case, the spectrum is gapless. Following the Solyom's analysis ||, we can predict the 
most preferred type of instability occurring in the ground-state of the system, corresponding 
to the most divergent correlation function: charge density wave (CDW), spin density wave 
(SDW), singlet superconductivity (SS) or triplet superconductivity (TS). Before discussing 
the results, let us remark that the on-site and bond-site interaction constants have always 
opposite effects; they occur only in the combination 



or, in the half-filled case, through Yi or Y 2 with a similar structure. Consequently, a bond- 
site repulsion acts (in the upper band) as an effective attraction and its effect is enhanced 
by increasing the electron density n, a fact already used in the hole superconductivity 



IV. PHASE DIAGRAMS 



Y = U - 8* sin(7m/2) 



(33) 
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mechanism |H|. The value of Y (or 1*1,2) fixes the position, in the (W, V)-plane, of the 
intersection between the "spin-line" (g a = 0) and the "charge-line" (g = | ^3 1 ) - 

Let us first consider the half-filled case (n = 1/2). The obtained phase diagram has 
the structure presented in Fig. 3. The "spin-line" is vertical and the "charge-line" is 
horizontal; both of them delimit not only different phases, but also the strong coupling regime 
from the weak coupling one in each corresponding (charge or spin) sector . The critical 
value V2 = — Max {Yj., I2} /2 of the inter-site interaction constant separates a dominant 
superconductor region (V < V2) from a density fluctuation one (V > V2). Analogously, by 
decreasing the intra-bond interaction constant W below the critical value Wi = —Y/A, the 
triplet states (spin-density or superconductor) disappear. Depending on the values of Ua,b 
and Xa,b, the origin W = V = can be found in principle in any of the four regions, but 
the most probable case for a real system is that considered in the picture (Y > 0, Y 12 > 0). 

In the non-half-filled case, #3 is zero and consequently g p does not renormalize; thus the 
charge sector is always in the weak coupling regime and the "charge-line" now separates only 
different phases. Graphically, two things happen: (i) the intersection between the "charge- 
line" and the "spin-line" moves in such a way that the origin W = V = lies inside either 
the region II (for Y > 0) or the region IV (for Y < 0); (ii) by increasing the density n, the 
"spin-line" rotates anti-clockwise while the "charge-line" rotates clockwise (by decreasing 
n, the effect is reversed). Consequently, on starting from the half-filled case and increasing 
the electron concentration the regions I and III become smaller while the regions II and 
IV grow, as it is shown in Fig. 4 for Y > 0; the intersection of the "spin-line" with the 
coordinate axes are given respectively by Wi = —Y/4 and Vi = Y / [2 cos(irn)} while in the 
"charge-line" case by W2 = V/[4cos(7rn)] and V2 = — Y/[A + 2cos(7rn)]. The effect of the 
model parameters on the ground-state phase diagram can easily be inferred from the above 
discussion. 
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V. THE CASE OF CU0 3 CHAIN 



The use of the RG method implicitly assumes a weak coupling regime; the scaling equa- 
tions are obtained in the second order of the perturbation theory and consequently the 
coupling constants have to be small in comparison with the bandwidth. For the 3/4-filled 
undistorted CuO% chain we have t — t — 1.4 eV with a gap 2 A = 2e = 1.23 eV JTT] and 
from Eq. (5) it follows a bandwidth around 2.25 eV (and thus a quite large gap). On the 



other hand, Ua ~3 - 4 eV, Ub ~1 - 3 eV [ll|, i.e. we are in an intermediate coupling 
regime; in this respect, our estimations for the CuO% chain presented below could provide 
only a qualitative description of the system. The restriction to the partial filled (upper) 
band, justified also in the weak coupling limit, supports the same conclusion. 

Excepting the on-site interaction, the values of the other coupling constants for the CuO^, 
chain are less known. For a non-alternating chain, using a Kronig-Penny model for the ionic 
potential and a screened Coulomb inter-electronic potential, Campbell et al. [pT2] , p~3|| analyzed 
the relative magnitude of the Hubbard-type constants for various ratios of the overlap (of 
the Wannier functions) and screening parameters: in general U > V > X > W, but in the 
extreme screening limit U > \X\ > V — W. In most real problems, the effect of the exchange 
and pair hopping (W) can be neglected and we shall assume this fact as valid also for an 
alternating chain. The phase diagram corresponding to a half-filled upper band (see Fig. 
3) in this case (W = 0) and in terms of the original coupling constants (with A/B = 3/4, 
as for the CuO% chain) looks like in Fig. 5. The oblique heavy line from Fig. 5 is defined 
by Y\ = Y 2 [see Eq. (32)] (the dashed line corresponds to Y\ = —Y 2 ) and it is independent 
of the inter-site coupling constant V; its slope depends only on the ratio A/B. The two 
heavy halflines (horizontal and vertical) start from the point Q, sliding on the dashed line by 
changing V (upwards by decreasing V; at V = 0, Q coincides with P); they delimit the region 
where a superconductor state can occur. As can be remarked from Fig. 5 (drawn in the 
particular case Xa = Xb = —V = 0.2 eV) the region of SS and TS (denoted by I in Fig. 3) 
disappears for V > 0. The values of the inter-site and bond-site coupling constants for which 
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the 3/4- filled undistorted CuO$ chain is superconductor can be determined also from Fig. 
5; for example, in the particular case of V = 0, a superconductor state can be reached only 
if X A /U A > B/(A^2A) ~ 0.24 and X B /U B > A/(A^2B) ~ 0.13 [See Eq. (32)]. For a non- 
alternating chain, the same condition requires Y < 0, i.e. Xq/Uq > I /(Ay/2) ~ 0.18 |L4] [See 
Eq. (33), where we replaced the renormalized coupling constants U and X by respectively 
Uq and X of a non-alternating chain]. In order to get a superconductor state, we thus 
need a bigger ratio of the (A) bond- and on-site couplings in the alternating case than in 
the non-dimerized one; however, in average, the ratio is practically the same and we could 
conclude that, within the present formalism, the alternating structure does not play an 
essential role in the occurrence of superconductivity in this system. Let us also note that for 
a non-alternating structure and in the extreme screening limit (when V can be neglected), 
the ratio (computed from atomic orbitals) Xq/Uq has been estimated in the range 0.15 - 
0.18 fL3| ; and it is closed to the value required for a superconductor state. 

Finally, let us note that for realistic values of the parameters, there is a coexistence of 
CDW and SDW fluctuations in the system (see the shaded rectangle from Fig. 5). This 
phase could be interpreted as follows: any allowed density modulation has a wave vector 2kp 
and consequently for the half-filled upper band the length wave of the density fluctuation is 
equal to two lattice constants. From the alternating structure it follows we could imagine 
SDW and CDW as located on A and B sublattices respectively (a spin density wave requires 
a stronger on-site repulsion than a charge density one). 



VI. CONCLUSIONS 

The results of this paper can be summarized as follows: the canonical transformation 
diagonalizing the one-particle Hamiltonian for an alternating chain with two non-equivalent 
sites per unit cell has been used to find the expressions of the Hubbard-type interactions, 
initially introduced in the site-representation, in the corresponding two-band model; the 
consideration has been restricted to the particular case of only alternating on-site energies 
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(and equal hopping amplitudes) and a gap between the two bands large enough to take into 
account only the processes inside the partial filled band (and we chose the upper band, as 
for the 3/4-filled CuOz chain). A particularity of the alternating structure manifests itself in 
the dependence of the obtained potentials on the momentum conservation in a given process, 
i.e. if it is normal or umklapp one. The RG method (second order) has then been applied to 
this (one-band) model Hamiltonian; all the ^-constants have been obtained by evaluating the 
corresponding Hubbard-type interactions around the Fermi points. The resulting ground- 
state phase diagrams have been analyzed in terms of all coupling constants and for an 
arbitrary band filling. The effect of the bond-site interaction is to renormalize the on-site 
one. In the chosen coordinates (intra-bond, inter-site) and at half filling, a critical inter-site 
interaction controls the position of the "charge-line" (separating a dominant superconductor 
region from a density fluctuation one), and a critical intra-bond coupling determines the 
"spin-line" (separating a region of singlet states from one with possible triplet states); the 
effect of the density is to rotate these two lines and to change the position of their crossing 
determined by a certain relation between the on-site and bond-site coupling constants. The 
particular case of the 3/4-filled undistorted Cu0 3 chain has been discussed in a separate 
section, where the corresponding phase diagram has been drawn in terms of the original 
coupling constants (neglecting W term, usually very small); the values of the inter-site and 
bond-site couplings required to get a superconductor state (according to our approximate 
description) have been also determined. 
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TABLES 

TABLE I. The values of the g-constants for the upper band of a chain with alternating on-site 
energies and Hubbard- type interactions (last column), compared with the similar quantities for a 
non-alternating structure (middle column). U, V, X and W are defined in the text by Eq. (21) 
and they are analogous, for our model, of the Hubbard-like coupling constants Uq, Vq, Xq and 
Wo from the usual case. By n we denoted the electron filling of the band, related to the Fermi 
momentum by akp = 7r(l — n) for the upper band of the alternating chain, or by akp = irn in the 
non-alternating case (0 < n < 1). 



non-alternating chain alternating chain 

g n 2V cos(27rn) + 2W -2V cos(vrn) + 2W 

g 1± U Q + 2V cos(27rn) + 8X cos(?rn) + 4W U - 2V cos(vrn) - 8X sin(vrn/2) + 4W 

g 2ll 2V + 2W cos(27rn) 2V - 2W cos(vrn) 

g 2± U + 2V + 8X cos(Trn) + 21L {1 + cos(2vrn)} U + 2V - 8X sin(7rn/2) + 2W{1 - cos(vrn)} 

g 3 (U -2V -4Wo)6 ntl/2 {{U A B A - AV2X A AB 3 ) - (A <- B)}5 n>l/2 

g 4 U + 2V + 8X cos(Trn) + 21L {1 + cos(2vrn)} U + 2V - 8X sin(7rn/2) + 2W{1 - cos(vrn)} 



TABLE II. The scaling of the considered model (14) to one of the two exactly soluble models - 
Tomonaga-Luttinger (TL) or Luther-Emery (LE), in each of the four regions defined by Eq. (27) 
in both the charge (c) and spin (s) sector. 





half filling (n = 1/2) 
c 


s 


not half filling 

c 


' (n + 1/2) 

s 


I 


TL 


TL 


TL 


TL 


II 


LE 


TL 


TL 


TL 


III 


LE 


LE 


TL 


LE 


IV 


TL 


LE 


TL 


LE 
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FIGURES 



FIG. 1. Alternating chain with two non-equivalent sites per unit cell. 

FIG. 2. The two bands corresponding to the one-particle tight binding Hamiltonian for an 
alternating chain with two non-equivalent sites per unit cell. 

FIG. 3. The phase diagram corresponding to a chain with alternating on-site atomic 
energies and Hubbard-type interactions for a half-filled upper band. Wi = — Y/4 and 
V2 = — Max {Yi, Y2} /2 [see Eqs. (32) and (33) in the text]. The response functions corresponding 
to the phases shown in the parentheses have a lower degree of divergence than the others. 

FIG. 4. The phase diagram of the considered model away from half filling. W\ = —Y/4 
and Vi = Yj [2 cos (7m)] determine the "spin-line", while the "charge- line" is determined by 
W 2 = Y/[4cos(vrn)] and V 2 = -Y/[4 + 2cos(vrn)] (here Y > and the band filling n > 1/2). 

FIG. 5. The phase diagram corresponding to a half-filled upper band of an alter- 
nating chain in terms of the original couplings (site representation) and for a negligi- 
ble intra-bond term; P and Q have the coordinates A{2) 1 / 2 X A A/ B,A(2) 1 / 2 X B B / A and 
A{2) l / 2 X A A/B - AVA 2 /B 2 ,A(2) 1/2 X B B/A- AVB 2 /A 2 } respectively and tan0 = B 4 /A 4 . The 
picture has been drawn in the particular case ^4/5=3/4, as for the C11O3 chain where the on-site 
couplings are estimated in the shaded rectangle. 
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